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Abstract 

This paper concerns the kinetic limit of the Dirac equation with random electromag- 
netic field. We give a detailed mathematical analysis of the radiative transport limit 
for the phase space energy density of solutions to the Dirac equation. Our derivation is 
based on a martingale method and a perturbed test function expansion. This requires 
the electromagnetic field to be a space-time random field. 

The main mathematical tool in the derivation of the kinetic limit is the matrix- valued 
Wigner transform of the vector- valued Dirac solution. The major novelty compared to 
the scalar (Schrodinger) case is the proof of convergence of cross-modes to weakly in 
space and almost surely in probability. The propagating modes are shown to converge 
in an appropriate strong sense to their deterministic limit. 

1 Introduction 

The Dirac equation is the relativistic version of the Schrodinger equation and describes 
very fast electrons propagating in an electromagnetic field. In this paper, we consider 
the semiclassical limit of the Dirac equation when the electromagnetic field is random and 
time-dependent. 

The problem falls into the category of high frequency wave propagating in highly het- 
erogeneous media, which has been modeled by radiative transfer equation in many areas, 
e.g., quantum waves in semiconductors, electromagnetic waves in turbulent atmospheres 
and plasmas, underwater acoustic waves, elastic waves in the Earth's crust. Such kinetic 
models account for the multiple interactions of wave fields with the fluctuation of the un- 
derlying media. In the so-called weak-coupling limit we will consider, waves propagate over 
distances that are large compared to the typical wavelength in the system and the fluctua- 
tions have weak amplitude with correlation length comparable to the wavelength. Most of 
the derivations of radiative transfer equation are based on formal expansions, e.g., a sys- 
tematic method to derive kinetic equations from symmetric first-order hyperbolic systems, 

*Department of Applied Physics & Applied Mathematics, Columbia University, New York, NY 10027 
(yg2254@columbia.edu; gb2030@columbia.edu) 



1 



including systems of acoustic and elastic equations, in the weak-coupling limit has been 
presented in [22] and extended in various forms in [1, 3, 15, 21]. 

Mathematically rigorous derivations are notoriously difficult in the setting of spatially 
varying randomness. Most proofs have been obtained for the Schrodinger equation with 
time-independent Gaussian potential, see [11, 23], with an extension to (discrete) wave 
equations in [17], and are based on the Neumann series expansion for the solution to the 
Schrodinger equation and appropriate estimates that allow passage to the limit. Such tech- 
niques were also used in [5] to obtain the random mixing of the phase of the Schrodinger 
solution and in [2, 24, 25, 16] to analyze limits of solutions with low-frequency initial con- 
ditions and large, high-frequency, random potentials. 

The derivation of kinetic limits is much simplified in the setting of time-dependent ran- 
dom coefficients and does not require the aforementioned infinite Neumann series expansion. 
Assuming a Markovian structure for the random potential enables us to use a martingale 
method and a suitable perturbed test functions expansion. A limit theorem for one di- 
mensional waves where such methods are used is given in [18] and more general ones in 
[9]. The same approach has been applied to Schrodinger equation in different settings, see 
[6, 7, 12, 14, 19]. 

The markovianity of the random coefficient simplifies derivations but is not necessary 
as was shown in [12]. In [10], a geometric approach is applied to more general initial data 
than in [11]. A renewal of the random field is used to get the appropriate mixing properties 
during the evolution. In [20], the authors consider random potentials that are correlated in 
time with finite range. They proved the semiclassical limit of Schrodinger equation by pure 
PDE techniques. 

For systems of equations such as the Dirac system or the linear hyperbolic systems con- 
sidered in [22], no rigorous results have been established even for time-dependent potentials; 
see also [4] for a recent review on the derivation of limiting models for wave and particle 
propagation in random media. In this setting, the energy matrix has to be decomposed into 
different propagating and non-propagating modes, which behave quite differently in the 
high-frequency limit. We focus here on the Dirac equation although most of the derivations 
carried out in this paper extend to the framework of linear hyperbolic systems, which will 
be considered elsewhere. It should also be mentioned that the energy bands in the Dirac 
equation are degenerate, which is distinct from the scalar (Schrodinger) case. 

The rest of the paper is organized as follows. We describe the problem setting and 
state our main results in Section 2. Next, we sketch the outline of the proof and discuss 
in detail the construction of test functions in Section 3. In Section 4, we then prove the 
convergence of different modes. The case of slower fluctuations in time, which allow us 
to recover limiting kinetic models with elastic scattering, is briefly discussed in Section 5. 
Conclusions and some further discussions are finally presented in Section 6. 

2 Main results 

In this section, we first describe our setting, next introduce the so-called Wigner transform 
as our main tool and then present our construction of the random electromagnetic field. 
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We finally state the main results of the paper. 

The Dirac equation in three dimensions of space reads 

ed t V £ + P(x, eD x )V E = 0, i£8 3 ,teR, (2.1) 

with the differential operator P(x, D) defined as: 

P(x, 0=i(j2 Wfo " eA k (x)) + m c 7 - eA (x)I^j . (2.2) 

Here ^ £ = ^ £ (t,x) £ C 4 is the wave function, e = h is our small parameter and stands for 
the Planck constant, e is the unit charge, mo is the electron's rest mass, c is the velocity 
of light, and A^x) Gl,fe = 0,l,2,3 are the components of the prescribed electromagnetic 
field. In particular, Aq is the electric potential and (Ai, A2, A%) is the magnetic potential 
vector. 

7 fc G C 4x4 , k = 0,1,2,3 are the 4x4 Dirac matrices, which are closely related to the 
2x2 Pauli matrices. Their elements are 0, 1, i, and they satisfy 

7 o* = 7 o ; y.* = _ ry k tk = 1} 2> 3; ( yy } * = 7 o 7fc) 

7 m 7 n + 7 n 7 m = 0, m ± n, 7 ° 2 = 7 4) 7^ = ~h, k = 1,2, 3. 
In our proofs, the explicit form of the Dirac matrices does not play any special role. 

The relativistic current density J e is a 4 dimensional vector with elements JjF given by 

Jl = ^*/ 7 fc ^, 

and the relativistic position density n(t, x) is given by 

n £ {t,x) = 4(t,x) = y £ *(t,x)^ £ {t,x). (2.3) 

We investigate the limiting behavior of the solution to (2.1) as e — > 0. Given the 
conservation of f Md n £ (t,x)dx and the fact that n £ (t,x) does not admit a closed-formed 
equation, the quantity we are interested in proves to be the Wigner transform of ^ /£ , which 
we introduce next. For more results on the Wigner transform, see [13]. 

2.1 Wigner transform and pseudo-differential calculus 

We introduce the Wigner transform and some pseudo-differential calculus that is to be used. 

The matrix-valued Wigner transform of two spatially-dependent d dimensional vector 
fields u(x) and v(x) is defined as 

W[u,v](x,0= [ j^e^u&x-^yix + ^dy, (2.4) 
J R d (2ir) a 2 2 

where v* is the transposition and possible complex conjugate of v. It may be seen as the 
inverse Fourier transform of the two point correlation function of u(x) and v(x), where we 
define Fourier transform using the convention /(£) = f Rd e~ l ^' x f(x)dx. We check that: 

(W[u,u](x,0T = W[u,u](x,0- (2.5) 
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We also verify that 

W[u,v]{t,x,Z)d£ = (uv*)(t,x), (2.6) 



and this allows us to interpret the Wigner transform as the energy density in phase space, 
although the Wigner transform is positive only in the limit e — > [13]. 

We recall from [1] some simple results about pseudo-differential calculus that are needed 
in this paper. 

Proposition 2.1. let P{eD) be a matrix-valued pseudo- differential operator defined by 

P{eD)u(x) = [ e^Ptiep)^)*. 

Then we have 



W[P(eD)u,v}(x,C) = P(i£ + ^)W[u,v](x,0, 
and 

W[u, P(eD)v] (x, e) = W[u, v] (x, £)P* (t£ - 

where W[u,v](x,£)P* (i£ — sD/2) is defined as the inverse Fourier transform Fp\ x of the 
matrix W[u,v](p,£)P*(i£ — iep/2). 

Proposition 2.2. let V(x) be a real matrix-valued function, then we have 



W[V&u,v](x,0= / e -— [ V(p)W[u,v](x,C-^)dp, 



and 

a ix-p/e 



W[u,V(-)v](x,0= / e -— I W[u,v]{x^ + ^)V\p)dp, 
e J R d (27r) a 2 

where V is the Fourier transform of V component by component. 

Remark 2.3. Throughout the paper, we assume that spatial dimension d = 3 although we 
occasionally use "d" for expressions that hold independent of dimension. 



2.2 The random field 



The Ak{x) appearing in (2.1) are components of a random electromagnetic field, which we 
assume to be time-dependent and have mean zero. The non mean zero case can be handled 
similarly. 

We follow the same construction of the random field as in [6, 7). V is the set of measures 
of bounded total variation with support inside a ball Bi = \p\ < M: 

V = {V = (Vb, Vi, V 2 , V 3 ) : f \d%\ < C, supp V k C B L , V k {p) = V k *(- P ),k = 0, 1, 2, 3}. 
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Let A(t,p) = (Ao(t,p)dp, Ai{t,p)dp, A2(t,p)dp, As{t,p)dp) be a mean-zero Markov process 
on V with generator Q. The time-dependent random field Ak(t,x) is given by 

and is real and uniformly bounded. We assume Ak(t,x) is stationary in t and x and the 
correlation functions are defined by 

R mn (t, x) := E{A m (t + s, x + y)A n (s, y)},m, n = 0, 1, 2, 3. 

Furthermore, we have 

E{A m (t + s,p)dpA n (s,q)dq} = (2n) d Rm n (t,p)6(p + q)dpdq, 

where the power spectrum R m n is the Fourier transform of R mn (t,x) in x: 

R mn (t,p)= R mn (t,x)e- ix - p dx. 

For simplicity we assume that R mn (t,p) e S(R x R d ) and define the space-time Fourier 
transform R mn {ui,p) as 

R mn (u,p) = / R mn (t,p)e- iujt dt. 



We assume that the generator Q is a bounded operator on L°°(V) with a unique invariant 
measure tt(V), i.e., Q*7r = 0, and there exists a > such that if (g,7r) = 0, then 

lle'SllL-(V) < C\\g\\ Loo{v) e~ at . (2.8) 

The simplest example of a generator with gap in the spectrum and invariant measure tt is 
a jump process on V where 

Q9(V) = [ g(V x )dK{Vi) -g(V), [ dv{V) = 1. 

Given (2.8), the Fredholm alternative holds for the Poisson equation 

Qf = 9 

provided that (g,ir) = 0. It has a unique solution / with (/, vr) = and ||/||l°°(v) < 
C||5 , IIl oo (V)- The solution is given explicitly by 

oo 

rQ, 



f(V) = - / dre rU g(V), 
Jo 

and the integral converges absolutely because of (2.8). 
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2.3 Main theorem 



Before stating the main results we first derive the equation satisfied by the Wigner trans- 
form. 

Recall the equation 

ed t V £ + P(x, eD x )V £ = 0. 

We replace Af.(t,x) with y/eAk(-, -) in the weak coupling limit, so the equation for ty 6 
becomes 

ed t ^ + i (£ 7 V(^fc - V~eeA k (-, -)) + m c 7 - V~eeA (-, -)lA * £ = 0. (2.9) 

Define W e (t,x,£) = W[V £ (t, .),* £ (i, .)](»,£)> and b y Proposition 2.1, 2.2 we can derive the 
equation satisfied by W £ 

3 T— . 3 T— V 

£6^+ ]T c 7 V^(^ + ^-)W 6 + W B P fc *(i£ " ^-)c 7 V + im c 2 (7 ^ £ - W el °) 

k=l k=l 

3 3 

fc=0 fc=0 

(2.10) 

where the symbol Pk(£,) = £fc; anci the operators JC £ ,K £ are defined as 

X*A(t 1 x,0= / ^%^e^ £ A(^,e-f), 
J M d (2ir) a 2 

J R d (27Tj fl 2 

We define the L 2 inner product for matrix valued functions as follows: 

(F(y),G(y)) = j Tr(F*G)(y)dy. (2.11) 

The Dirac equation (2.9) preserves the L 2 norm of ^> £ , and by the property of Fourier 
transform, we have Vi E [0, T], 

||W £ (t)|| L2(RM) = (2Trer d/2 \\* £ (t)\\ 2 LHRd) = (27r e )- d / 2 ||^(0)||i 2(R(i) . (2.12) 

Thus by choosing appropriate initial condition \& e (0) such that 11^(0)11^2^) is of order e d / 2 , 
we have ||W / £ (t)|| L 2( K 2d) uniformly bounded for t £ [0, T] and A £ V. Moreover, we assume 
that We(0, x, £) converges weakly as e — > to Wo(Q,x,£) € L 2 (M. 2d ) without restriction to 
a subsequence. 

Expanding W £ = Wo + \fsW\ )£ + eW2, £ + ■ ■ ■, dividing (2.10) by e and expanding in e, 
by setting the 1/e order term to be zero, we obtain that 

3 3 

ci°i k PM)Wv + WoY, Pk(iOH°l k + imo^i^Wo - W ol °) = 0. 
fc=i fc=i 
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This enables us to define the dispersion matrix 

3 



(2.13) 



k=l 



and we expect the limit to satisfy that 

Q(OW (t, x, = W (t, x, 0Q(0- (2-14) 
By the property of 7°7 fc , we know Q is Hermitian, and Q 2 = (m^c 2 + |£| 2 )i4. Let 

\ ± (0 = ±^m 2 c 2 + \Z\ 2 , (2.15) 

be the eigenvalues of Q corresponding to the energy levels of electrons and positrons, re- 
spectively. The orthonormal eigenvectors are 







( 



yi (0 



V2A+(A + -£ 3 ) 
V2A + (A + -g 3 ") 
2A + 

mpc 



aT^ 

2A + 

V V 2A +( A +"«3) / 



2A+ 



2/2 (0 





V2A + (A + ^ ) 
2A+ 



rape 

\ V 2A +( A +-«3) / 



where x\,X2 correspond to A + and 2/1,2/2 correspond to A_. 

Remark 2.4. The explicit form of the eigenvectors do not affect our results. 

Because {x±, X2, yi, 2/2} is a basis in C 4 , we can decompose W £ as 



(2.16) 



On one hand, we have W* = W £ , so 



1,2 a 



12 



*21' 



1,2 612 = 6^. 



On the other hand, since W e is uniformly bounded in L 2 , we have that all the coefficients 
a?-, c?-, dlj belong to C([0, T]; L 2 (R M )). The dispersion relation (2.14) implies that c?-, d\ 
should converge to zero as e — > 0. Thus we write the limit 



«J=1,2 i,i=l,2 



(2.17) 
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We denote by II± the orthogonal projection of C 4 on the eigenspace associated to A±, and 
they are 

11+ = xix\ + x 2 x* 2 = \{h + -p-), 
z A-|- 

II- = y\y\ + 2/22/2 = \{h - -^-)- 
Define a% = Ti(U + W E ) = of x + a| 2 , ai = Tr(n_W £ ) = 6f x + 6 22 . We check that 

(ai + ai)(f,x,£)d£ = n £ (t,x). 



/ 



Therefore we can interpret a+ and ai as energy densities in the phase space corresponding 
to electrons and positrons, respectively. They are the physical quantities we are interested 
in. 

Define L 2 M (R 2d ) = {f\\\f\\ L 2 {R 2 d) < M}. The space we use for (cf + (t),cf_(t)) is 

H M = iM(» M )$4(R M ) = {X = (X 1 ,X 2 )\X 1) X 2 G L 2 M (R 2d )} 
endowed with the metric 

wr) .^ W-Ml t K-MI \ 

n.= l ^ ' 

where {e n } is orthonormal basis of L^IR 2 ^). It is straightforward to check that dn induces 
the weak topology on L 2 I (M 2d ) L 2 I (M 2d ) and (Hm,cIh) is a complete separable metric 
space. We have the following tightness criteria about process taking value in -Hju- 

Proposition 2.5. Suppose P e is the family of probability measures induced by (ai,ai) on 
C([0,T};H M ), if for any f G C 1 ([0,T];S(R 2d )), the process (f,a e ±) G C([0,T];M) is tight, 
then P £ is tight. 

Now we can state our main theorem. 

Theorem 2.6. Suppose W £ solves (2.10) with an initial condition W e (0,x,£) that is uni- 
formly bounded inL 2 {M? d ) and converges weakly inL 2 (M? d ) to a fixed, deterministic function 
W (0,x,£) G L 2 (R 2d ) as e ->■ 0. We decompose W £ as 

W e= Y "■r r ' r ', ■ Y h U^'"'i ■ Y <$j x iVj+ Y '' ■!!>'■'' J' 

»J=1»2 «J=1,2 i,j=l,2 ij=l,2 

T/ien /or i/ie cross modes cfj , dfj we have 



V (\ / I &(t,x,£)f(t,x,Odxd£dt\ + \ f [ dUt,x,0f(t,x,0dxd^dt\] <C /iTX /i 

(2.18) 
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for any function f £ C 1 ([0, T]; S(M. 2d )) almost surely. For the propagating modes, define 
a e + = af 1 + a|2, ot e _ = b\ l + 6§2 an d suppose P e is the family of probability measures induced 
by (afLaf.) on C([0, T]; Hm). Then, as e — > 0, ¥ £ converges weakly to the probability 
measure P = dr a + , a _), where (a+,a_) is £ae unique deterministic solution to the following 
transport equation system 

c£ ■ V x a + 

/t}° (2-19) 
A (|) = 

Tae initial conditions are given by a±(0,x,£) = Tr(II±Wo(0, x, 4)), and t/ie scattering op- 
erator T is defined as 

-,2 



3 

/ («-(<?) - «+(£)) y^ fe (£,g)i? fefc (cA + (g) + cA + (£),g- 4)dg, 



e 

+ 



(2vr)< . 

where we have g) + a)fc(4,g) = 1 and 

, t \ A + (g)A + (e)+69i+6^2 + 693 + ^c 2 

wi(4,g) 
^2(4,9) 
^3(4,5) 







2A + (g)A+(4) 


A, 




(4) + 4igi - 42^2 - 43^3 - mlc 2 






2A + (g)A + (4) 


A, 




(4) - 4igi + 42^2 - 43?3 - mlc 2 






2A + (g)A+(4) 


A. 


h(?)A- 


h(4) - 4igi - 4202 + 4353 - m\c 2 



2A + (g)A+(4) 



Remark 2.7. We have assumed that R mn = when m 7^ n. Similar results could be obtained 
for the general case with more complicated expressions we do not reproduce here. 

Remark 2.8. We see from the structure of the scattering operator that due to the temporal 
regularization, the energy ^fra^c 2 + |4| 2 is no longer conserved. Scattering is inelastic and 
we observe a coupling between the propagating modes a+ and a_. 

Remark 2.9. The uniqueness of solutions to the transport equation system comes from the 
fact that dt((a + , a + ) + < under the dynamics of (2.19). To see this, we only 

have to note that 

c£-V x a± 



9 



and 

(a+,T(a+,a_)) + (a_,T(a_,a + )) 



2 /• 3 

/ [(a+(g) - a+(£)) 2 + (a_(g) - a_(0) 2 ] q)R kk (c\ + (q) - c\+(0,q- 0<Lqd£dx 

3 

/ (a+(£) - a-(?)) 2 VwtK, q)R kk (cX+(q) + cA+(f )> 9 ~ Qdqd£dx < 0. 



e 2 



(27T) 



2.4 Comments on the limiting equation 

Using the same approach, we can more generally show that the matrix-valued process 
(A £ (t), B e (i)) has the weak limit (A(t),B(t)) satisfying a transport equation system, where 
we have defined A £ (t) = (a?.(t)),B £ (t) = (&?•(*)) and A(t) = ( aij (t)),B(t) = (6y(t)). In 
other words, we have a limiting transport equation system for af^, b\y 

Using the operators defined in (3.1), we formally have 

8 t W = A*W , (2.20) 



where Wo = Yli,j=i2 a ij x i x j + Yl% j=i 2 ^ijUiUj ■ To derive the equation satisfied by 
or bij, we only have to write (2.20) as d t x*W Xj = x*A*W Xj or d t y*W yj = y*A*W yj 
and compute x*A*WoXj,y*A*Woyj respectively. It should be mentioned that Wo does not 
satisfy (2.20) b ecause the dynamics would inevitably generate modes of Xiyj^y^Xj. 

Therefore, we have the limiting coefficient matrices A(i) and B(t) satisfying the following 
transport equation system: 

q a if\ f x*iA*W oXl x* 1 A*W x 2 \ _ ( ylA*W m ytA*W y 2 \ 

° tA[t) ~ V x* 2 A*W QXl x* 2 A*W x 2 J ' ~ \ y* 2 A*W oyi y* 2 A*W y 2 J " 

So a + = an + 0-22 and a_ = b\\ + b 22 satisfy 

d t a + = x\A*Wqx\ + x 2 A*Wox 2 , 
d t a- = ylA*W oyi + y* 2 A*W y 2 . 

If we calculate these expressions explicitly, we recover (2.19) in Theorem 2.6. 

Prom (2.19) we see that the equation of a+ — an + a22 does not involve a>ij,bij when 
i 7^ j. This could be illustrated as follows. 
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First of all, we write the equations satisfied by an and 022 explicitly. 



(2vr) c 



d t au + 



A+(0 



E E ^(9)^(07°7 lb ^(?)^(?)7 7*xi(0^d? 

fc=0 i,j=l,2 

E E ^(?)^(07 7 fe 2/i(?)y*(«)7°7^i(0^^ 
fc=0i,j=l,2 

E«n^) [K(£)7°7 fc zi(<?)| 2 + l^(e)7 7^ 2 (c?)| 2 ] ^dg 

fc=0 

3 



E a ^) Ni(07Vj/i(9)| 2 + |xJ(07Vy2(?)| 5 



l# fc dg, 



fc=0 



(2.22) 



(27T) C 



<9 t a 2 2 + 



eg • V x a 2 2 



+ 



fc=0i,j=l,2 

E E ^■('?)^(e)7 7 fc y J (9)y;(9)7 7^ 2 (e)^^ 

jfc=0i,j=l,2 
3 



fc=0 

E°22(0 



^(07°7^i(g)r + |^(e)7 u 7 fc x 2 (g)| 



R kk d( l 



xmi Q i k yM\ z + l^(07 u 7 ft y 2 (<?)| 



fc=0 



(2.23) 

where i?^. := R kk {c\+{q) - c\+(£),q - 0, R tk := Rkk{c\+(q) + cA+(g),g - £). So the 
equations for an and 022 are indeed related to a\2, 0,21, 612, &2i- However, if we sum (2.22) 
and (2.23), with the observation that 

7VQ(07°7 fc Q(<?) + QtahV QtehV = c k h 

for some constant c&, we have all terms concerning ay, 6y with i 7^ j cancel out. 
Take ai2 for example, we note that the coefficient of ai2 is 

(I) =x* l (0l°l k x l (q)x* 2 (q) 1 1 k x 1 (0 + ^(OtVM^MtV^CO 

=x* 2 ( q ) 1 ° 1 k (x l (Oxi(0 + x 2 (Oxmh i k Mq) 

x 2 (q)l°l k \{h + #^)7°7^l(9) = x* 2 (q)l(h + tV^tVM?) 



^2(9)Q(9)7°7 fc Q(e)7°7^i(g)/(2A + (g)A + (0) 
^(g) 7 7 fc Q(67°7 fc Q(9)xi(g)/(2A + (g)A + (0), 



A+(0 



so 



2(1) = [ 7 °7 fc Q(07°7 fc Q(Q) + Q(?)tVQ(07VJ xi(?) = 0. 

Similar expressions hold for a 2 i,bi 2 ,b 2 \. 
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3 Outline of the proof 



The proof of cross modes converging to zero comes from the equation itself and the a 
priori L 2 boundedness of W £ . For the propagating modes a§_, the basic idea is to use the 
Markovian property of A^(t,p)dp, k = 0,1,2,3. We construct appropriate test functions, 
then prove the approximating martingale property which will be specified later. Using the 
approximating martingale inequalities, we can prove the uniqueness of the limit; together 
with the tightness of (afi, a £ _), we finish the proof of the main theorem. 

First, we define the operator C and A 

.2 3 



Z u poo p 

^0(^6=7^ E / / 7% m ^n(r,g)e-^«-2^)[Ao^-g)7V-7 7 n Ao(6]e- rA2(e ---«)^dr 

W m,n=0 Jo J ® d 
2 3 poo p 

~7Tu E / / e- J - Al «-i^)[A (67 7 m -7 7 m Ao(e-g)]e-^^-f^ 7 V^(r J g)^r, 

1 3 1 

^Ao = - E( c 7°7 fc A-Ao + -I^AoctV) + £A , (3.1) 



where 



k=l 

3 

k=l 

3 

Mt,P) = E P <^ - f )c7°7 fc +^ c 2 7 c 
k=l 



(3.2) 



We will see later that for any test function Ao such that QXo = XoQ, we have (Ao, Wo) 
formally satisfies the equation 

d t (\o,W ) = ((d t + A)X ,W ), 

and this motivates us to define the first approximating martingale functional 
G{ o : C([0,r];L 2 (R 2d )) C([0,T]) as 

Ol [W](t) = (A ,W)(i) - I ds((d t + A)Xo,W)(s). (3.3) 

J o 

Assume P e is the probability measure induced by W £ on C([0, T\ \ L 2 (M. 2d )). Our first goal is 
to show that under P e , Q\ [W]{t) is an approximating martingale. More precisely, we will 
show 

\& % {G\ [W]{t)\F s } - g\ [W]{s)\ < C Xo , T V~e, (3.4) 
uniformly for all W G C([0, T]; L 2 (M. 2d )) and < s < t < T. 

Supposing (a+,af_) =>• (a+,at-) and choosing appropriate test functions Ao in (3.4), we 
obtain from the above bound that (E{a+}, E{a_}) satisfies the transport equation system 
(2.19) by setting s = and e — > 0. 
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On the other hand, for any test function Fq and Go satisfying the dispersion relation, 
i.e., QFq = FqQ,QGo = GoQ, let fio = Fq (g) Go and define the second approximating 
martingale functional Q 2 ^ : C([0,T]; L 2 (R 2d )) ->■ C([0,T\) as 

Gl [W](t) = (F ®G ,W®W){t)- [ ds(F ®[{dt + A)Go] + [{d t + A)F ]®G ,W®W)(s). 

Jo 

(3.5) 

Our second goal is to show that under P e , ^ [W](f) is an approximating martingale, i.e., 

\^{G 2 jW](t)\F s } - G 2 jW](s)\ < C w ,tA (3-6) 

uniformly for all W £ C([0, T]; L 2 (M 2d )) and < s < t < T. 

Similarly, by choosing s = and letting e — >• 0, we can show that E{a+ £5 a+} = 
E{a+} ® E{a + } and E{a_ ® a_} = E{a_} E{a„}, so (a+,a_) is deterministic hence 
the solution to the transport equation system (2.19). 

We next go through the construction of test functions. 
3.1 Construction of Ai e , A2, e 

In order to obtain the approximating martingale inequalities, we have to consider the con- 
ditional expectation of functions F(V,W). The only functions we are interested in are 
those of the form F(V, W) = (X(V),W). Given a function F(V, W), if we denote by P e 
the probability measure induced by (W £ (t),A(t/e)) on the space C([0, T]; L 2 (R 2d )) x V, the 
conditional expectation is defined as 

^wy tt {nV,W)}{r) = E^{F(A(r),W(r))\A(t) = V,W(t) = W},r> t. 
The weak form of the infinitesimal generator of the Markov process induced by (W e (t),A(t/e)) 

is 

J-E^{<A(t>), W)}(t + h)\ h=0 = (- £ QX + d t X + A £ X, W), (3.7) 

provided we write the equation (2.10) as dfW £ = A^W £ with A* the adjoint of A £ . 
We have 

AeX = --j2 + - - E \Pk{H - ^)c 7 °7 fc + -*m c 2 7 °A " -^ c 2 A 7 ° 

& Zt & Zt s & 

k=\ k=l 

~ ~U E tV*2A + -U E ^A 7 °7 fc . 
v k=o * k=0 

(3.8) 

Let A e = A + VeAi i£ + eX 2 , £ with A satisfying Q(£)Ao(£) = A (£)<2(0- Plugging this 
expressions into (|Q + dt + Ar)A e and equating like powers of e, we find that the term of 
order 1/e equal to zero. 
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Considering the term of order 1/y/e, we introduce the fast variable z = x/e and define 
Ai, e = X\(A, t, x, x/e, £), where Ai = X±(A, t, x, z,£) solves 

QXi - A 1 X 1 - X\A 2 = F u (3.9) 

with Ai = T z -*p\i, A\,Ai defined in (3.2) and 

F 1 = ieG t (p)\ (£ - |) - ie\ (d + |)G t (p), 
3 t 

G t ( P ) = Y J i°i k M-,p)- 

k=0 

The solution to (3.9) is given by 



e 



iz-p roo 



___ I e rQ e- r ^F ie - rM drdp. 







Similarly, we can define A 2>e = X 2 (A,t,x,x/e,^) with A2 = \2(A,t,x,z,£) solving 

QA 2 - A X X 2 - X 2 A 2 = F 2 + £X {2ir) d 5(p), (3.11) 

where 

The solvability comes from the fact that E{i^} + £Xo(2ir) d 5(p) = 0, and the solution is 
given by 

A 2 = -/ TTTd / e rQ e- rA i(F 2 + £Xo(2n) d 5(p))e- rA *drdp. (3.12) 

By (3.7), we know that 

Gi £ [W](t) := (A e , W0(t) - ! ds((-Q + d t + A e )X e ,W)(s) 

Jo £ 

is a P e -martingale. With \ e = Ao + y/eXi f£ + ^A2 )£ , we have 

Q%m® = oi [w](t) + (c u w)(t) + f 'ds(c 2 + c 3 + c 4 ,w)(s), 

Jo 

where the correctors Cj are 

C\ = \/iAi )E + eA 2j£ , 

C 2 = -cMVeAi e + eA 2j£ ), 



C7 3 = ^e^7°7 fc ^A 2 , £ - v^e £ ^A 2 , e7 V, 

fc=0 fc=0 

3 3 



k=l k=l 
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Remark 3.1. In C4, the derivative D x is with respect to the slow variable. 

To prove that C/J[ [W](£) is a ^-approximating martingale, we need to show that the 
following correctors are small: 



i=i i=i 



uniformly in t G [0, T]. 



3.2 Construction of /xi e , ^t2,e 

First of all, we derive the equation satisfied by W £ (8> W £ (t, x\, x 2 , £1, £ 2 )- Define the 16 x 16 
matrices 

r. \ / 7 * \ 



We have 

3 



/ 71V4 


7l2*4 


721*4 


722-^4 


7*1*4 


732*4 


V 74*1*4 


742^4 


3 r 





rl 3 h 7 3 fe 4 

.,k t -,k 



7 fc 
7 fc 
\ 7 fc J 



k=l 



W £ ®W £ + im c 2 (r" + T")W £ ® w £ 



fe=l 



cr°r*jfc(« - ^) + cf f fc p fc * 2 (;e - ^) 



imotHWe® w e (r u + r°) 



-ie^^(r°r fc /C^ + f°f k IC k 2 )W £ ®W £ + ie^W £ ® W £ ^(T ^*! + f°f fc /C^ 2 ) = 0, 

fc=0 fc=0 

(3.13) 

where P k i, P^, JC k { , JC k { are the corresponding operators with respect to (xi,£i),i = 1,2. 
The weak form of the infinitesimal generator of the Markov process (/x, W £ (g) W e ) is 

J-Ej^MVO, W (8 W0}(« + h)\ h=0 = {{- £ Q + d t + B e )fi, W®W) 
if we write (3.13) as d t (W £ <g> W £ ) = B*(W £ ® W £ ). Define /x £ = /i + v^i.s + e/x 2)e . Then 

G'ftWW) ■= <Afc, W ® " T + & + ^ ® W )( s ) 

Jo e 

is a P e -martingale. 

Consider (-Q + (9 t + B £ )fi £ , after expanding in e, the term of order 1/e is 

( J ) = " i E c [ r ° r "^i«) + f °f fc p,: 2 (ie)] /10 - \i*> E c [r°r fc P fc i(^) + f °r k p k2 ^) 



fc=i 

+ ^m c 2 (r° + f °)/i - ^im c 2 /Uo(r° + f°). 



fc=i 
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If we choose fj,o = Fq ® Go with Fq, Gq both satisfying the dispersion relation, i.e., 

Q(OFo(0 = F (OQ(0,Q(OGo(0 = Go(OQ(f). 

we can check that (I) = 0. For the term of order 1/y/e, similarly we introduce fast variables 
z\ = xi/e,Z2 = X2/S, and define 

~ ~ X\ X 2 

/ii j£ (^,i,xi,x 2 ,6,6) = fJ-i(A,t,x 1 , —,%2, 



with Hi solving 

3 



fc=i 



/zi + im c 2 (r° + f °)/zi 



fc=i 
3 



r°r fc p fel (^-^) + f°f fc p fe2 (^-^ 



im cVi(r° + f ) (3.14) 



=ie ^(r°r fc /c^ + f °f k ic k 2 )no - ieno ^(r°r fc £* 1 + f °f k ic k 2 ). 

k=0 k=0 

We can check that the solution \i\ = F\ <g) Gq + Fq ® Gi, where F\,G\ solve (3.9) with Ao 
replaced by Fq,Go respectively. 

Remark 3.2. In (3.14), we replace xi/e by z% in the operators fc k ei ,K k ei . 
In the same way, we define 

f J L 2>£ (A,t,X 1 ,X 2 ,£ 1 ,&) = H2(A,t,X!, —,%2, ~ 



and fj,2 solves 
3 



k=l 



T°T k PZM + ^-) + f°f fe P fc * 2 (i£ + ^) 



k=l 

3 

k=0 



r r k p kl (iz-^) + r f k p k2 (it-^) 



fi 2 + im c 2 (T° + f°)/i2 
im c 2 fi 2 (T° + f°) 



(3.15) 



^(r°r fc /c^ + f°f fc /c^ 2 )^i - iem J^iv^ic^ + f °f fc ^ 2 ) + l m . 



k=0 



We define L/xo such that the RHS of (3.15) has mean zero so the equation is solvable, i.e., 

3 3 

L/i = E{-ie^(r°r fc /c^ + f°f fc /c^ 2 )/ii + ie^i ^(r ^^ + f °f fc ^ 2 )}. 

k=0 k=0 

We decompose Ii/iq into two parts, i.e., hfi = L1//0 + L 2 ^o> where 

3 33 3 

Li/io =*eE{- r r fc /C^F! ® G + F x G ^ r°r fc ^ - f°f fc /C^F <g> Gi + F ® Gi ^ f°f fc ^ 2 } 



fc=0 



fc=0 



k=0 



k=0 



--(£F )®G + F ®(£Go), 
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and 



L 2 /i = ieE{- J2 r^J&FoSGi+FoSGi r !^-^ f °f ^F^Gq+F^Go £ f °f fe /C £ fe 2 }. 



fc=0 



fc=0 



fc=0 



fc=0 



By this decomposition, we can write // 2 = Fo ® G 2 + F 2 (g) Go + jj>2 > where F 2 , G 2 solve (3.11) 
with Aq replaced by Fo, Gq and /2 2 solves the following equation 



Q#2 - X] C 



fc=l 



r°r fc p^(ie + ^ ) + f °f k p* k M + ^ 



k=l 



/i 2 + zmoc 2 (r + f°)/x 2 

im c 2 /i 2 (r + f ) 

3 3 

=L 2Mo + ieE{J^ r°r fe X:^Fo ® Gi - F ® G x r r fc ^ + f °f k K k 2 F x ® Go - Fi <g> G ^ f°f A 

fc=0 fc=0 fc=0 

With fi e = fiQ + y/e^ii^ + e^ 2j£ , we have 

O^(*) = 0jo[ W ](*) + <^ W ® W O(*) + / ds(J2Ci,W®W)(s), 



k=0 



i=l 



and the correctors G, are 



Gi = \/e/ii )£ + e^ 2 , e , 

G 2 = -dt(y/£fJ,i,e + £fJQ,e), 

3 3 

g 3 = iev^^(r°r fc 4\ + f°f k K k £2 )fi 2 ,s - ieV~e^,e £(r r fe £ £ \ + f °f k t k 2 \ 

k=0 k=0 



fc=i 

+-(v / i/ii, e + £M2, £ )X] 



r°r fc F fc * 1 (^) + f°f fc F fc * 2 (^; 



fe=i 



r°r fc F fel (- £ f £ ) + f°f fe F fe2 (-^) 



Remark 3.3. In G5, -D^ is with respect to the slow variable. 

To prove that Q^ [W] (t) is a P £ -approximating martingale, we need to show that 

5 5 

i=i i=i 



uniformly in t G [0, T]. 
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4 Proof of the main theorem 



In this section, we first prove the convergence of the cross modes and then of the propagating 
modes. In the end, we prove the tightness result and finish the proof of the main theorem. 

From now on, we use the notation a < b when there exists a constant M such that 
a < Mb. 

4.1 Cross modes 

First of all, we show that the cross modes converge to zero weakly almost everywhere. For 
any test function F, we have 

dt(F,W e ) = ((d t +A £ )F,W £ ). 

Therefore, 

e(F,W £ )(t)-e(F,W £ )(0) = [ ds((ed t + eA £ )F,W £ ){s). 

Jo 

By (3.8), we rewrite eA £ F and have 

3 3 

eA £ F =(J2 c 7 V^fc + im c 2 7°)F - i£ fc c 7 V + *m e 2 7 °) 

k=l k=l 

3 3 3 3 

-e £ try VP** (§)F " eF £ P k (- f )c 7 V - v^e £ 7 V/C*F + yQe £ /C £ fc F 7 V. 

fc=l k=l k=0 k=0 

Since (W £ (i), W E (t)) < M uniformly in t £ [0,T] and A G V, /C^/C^ are bounded operators 
on L 2 (M. ), and Q(f ) = ^Li 7°7 fe £fc + ^oc 7 °, we have 

| [ T ds(F,QW £ -W £ Q)(s)\ = \ [ T ds(QF-FQ,W £ )(s)\<^. 
Jo Jo 

Recall that 

w £ = ^2 • ^2 h %yiy]+ c tj x iyj + ^2 d %y iX *v 

ij=l,2 M=i,2 »J=1,2 »,J'=1,2 

so QW e - W e Q = 2A+ Ei J= i, 2 tyiVj + 2X - E ij= i,2 If we choose F = fmyyK 

for some good function /, we have 

T ds(f,c! j )(s)\<^- 







The same is true for cf?- if we choose other test functions. Thus we conclude that the cross 
modes cf^dfj converge to zero weakly almost everywhere. 
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Remark 4.1. We do not necessarily have the weak convergence of cfj,dfj as processes in 
C([0, T]; L 2 (M. 2d )) here. The following is a heuristic argument. 

We rewrite the equation (2.10) satisfied by W £ : 



dtW £ + - 

e 



k=\ k=l 



+ ... 



1C 

=d t W £ + -{QW £ - W £ Q) + . . . = 0. 

(4.1) 

All the terms that do not show up here are those of order 1 or 1/ y/e. So we have 

x\dtWeVj + -x*(QWe - W £ Q) Vj + . . . = 0, 

which leads to 

2ic\ 



d t 4j + -^ c tj + • • • = o. 

We see that as e — > 0, the cross modes cfj are highly oscillatory in time. The same happens 
to dfj. Prom this perspective, we can only expect the weak convergence as a L 2 function 
rather than a process in C([0, T];L 2 (R 2d )). 

4.2 Convergence of the expectation 

We first define some notation used in the proof. For two matrix-valued functions A(y) and 
B(y), we have 

\A\ < \B\ ^V(m,n),\A mn (y)\ < \B mn (y)\ 

\A\ < max \B\ <3> V(m,n), \A mn {y)\ < max \B mn (y)\ 

\y\<M \y\<M 



\A\ < ||jB|| £ oo( V ) V(m,n), \A mn (y)\ < H-BmnWIlL^Cv) 

We have the following lemmas concerning Ai )£ , \2, £ - 
Lemma 4.2. (\i >e (t), \i, £ (t)) < C\ ,t uniformly for t G [0,T] and A G V. 

Proof. Recall from (3.10), we have \i iS (t,x,£) = \x(t,x,x/e,£), and 
\x(t,x,z,g) 

= - r dre* [ ^e-rMt,P) ie \ Gt (j>)\ (( - f ) - A (C + f )G»(p) 

Jo JR d ( 27r ) L 2 2 

Take one term, for example 

"oo r A^A.lt 



e -rA 2 (£,p)_ 



(I) 







dre rQ /" d PM^ e iz. Piee -rA^ P ) fc A ( g _ g )e -rA 2 (g, P )^ 
(2vr) d ' ' V ^ 2 J 
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Since Ak(-,p)dp is compactly supported and of bounded total variation, by (2.8) we have 

roo 

|(I)| < / dre~ ar max | e - rAl ^ p y 7 fc A (£ - ^) e - rA ^ p) \. 
Jo \p\<M 2 

Because Ai(£,p) = -icQ(£ + l),A 2 (£,p) = icQ(£ - §), any element in e ~ r ' Al ^ and 
e -rA 2 {i,p) j g uniformly bounded and thus 

| \ e -rM(^k Xo ^ _ P )e -r M ^ , < | [A (£ - f ) 

Since Ao is a good function, the proof is completed. □ 

Lemma 4.3. (A 2 , e (t), A 2 , e (t)) < C Ao>T /or t e [0,T] and IeV. 
Proof. The proof is similar with the one for Ai e . Recall that 

QA 2 - AiA 2 - X 2 A 2 = F 2 - E{F 2 }, 

where 



ic: 



Gt(q)M(p ~ q,£ ~ -)dq - j^— 
(2vr) a J R d 2 {27r) a 



Ai(p - <?,£ + -)G t (q)dq, 



and 



Ai(p,0 = - / dre rQ e- rAl( £' p *> 
Jo 



i e G t (p)Ao(e-|)-ieAoK + |)G t (p) 



-Ma(£,p) 



(27r )ir(#i --S 2 ), where 



therefore we can write F 2 - 

r roo 

B 1 = - dqGt(q) / dre rQ e- rAl ^-2' p - q) 
Jwi d Jo 



and 



G t (p - q)Mt ~ §) - A (£ + I - g)G*(p 



do / dre re e- rAl «+f G t (p - g)A (£ + «-£)- A (£ + ?)Gt(p - q) 



-rv4 2 (£-f,p-<?) 



-rA 2 ($+f,p- 



9) G t (g). 



Then we can write 



dre' 



dpe lz -P c _ rAl 



io V (2vr) d 

Considering the term B\ — EJBi, we have 



(2tt)< 



r(Bi -B 3 - £73!+ EB 2 )e- rAs . 



/* /'CO /• 

rQ / dpe lz -" p e- rM (B x -EB x )e- rM \< \ dre~ ar \\ dpe lz - p e~ rAl (B 1 - EB^e^ 2 \\ L o 
Jm d Jo Ju d 

roo r roo r 

< / dre- ar \\ / dpe iz -Pe- rM B l e- rM \\ L ^ {v) + / dre""! / dpe^e^EBie- 
Jo Ju d Jo il d 





oo 



'(V) 



li°°(V) " 
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|L°°(V)- 



Taking || J Rd dpe lzp e rAl B\e rA2 \\L°°(v) for example, we have 
/ dpe iz - p e- rAl ^ p) B ie - rA2 ^' p) 

JR d 

r r poo 

= - / dpe iz - p e- rAl ^' p) / dqGt(q) / dre rQ e- rAl ^-^ p ~ q) 
JR d JR d Jo 



G t (p - g)A (£ - |) - A (e + | - q)G t (p 



-rA 2 (£-&,p-q) e ~rA 2 {(,,p) _ 



Since G t (qr) = ££ =0 7 7^*(i, 9), we have 
| / dgl fcl (-,</) r dre rQ 

JR d £ JO 



e Jo 

t 



dpe iz - p e- rAl ^' p) -f -/ kl e- rAl ^-2' p -^ 1 1 k2 A k2 (-,p- g)A (£ - ^y- rA ^-\ ,P-q) e -rMZ,p)\ 

£ 2 

< f dq\A kl (-,q)\ r dre~ ar max \e^ A ^^e^ A ^-^h ^Hi ~ ^e^^^e-^^ 
jR d e Jo |p-«|<Af 2 

and so similarly we have 

| ^ e -rAi(^p) 7 7 fc le -rAi(C-§j»-g)^0 7 fc2^ o ^ _ P) e -rv4 2 (£-2 ,p-«) g-rA 2 (S,p) 



which completes the proof. □ 

By Lemma 4.2 and 4.3, we have ||(Ci, Ci)(t)||x,oo(v) < Cao.tE- For Cj,i = 2,3,4, the 
proof is similar so we omit the details here. In summary, we have 

4 4 

lK^a,^a)(t)ii LOO(v) <c Ao , Te 

8=1 8=1 

uniformly in t G [0, T], i.e., the 1? norm of Cj is of order y/e, so £/J[ [W](t) is an approxi- 
mating martingale under P e , i.e., 

\^ e {Gl [W](t)\T s } - Gl [W](s)\ < C Xo , T V~e. 

Now we can prove the following proposition about convergence of expectation. 

Proposition 4.4. If (a+,a £ _) converges weakly to some (a + ,a_) in C([0,T]; Hm), then 
(Ea+,Ea_) is the weak solution to (2.9). 

Proof. First of all, we point out that (Ea+,Ea_) has a fixed initial condition by our choice 
of W e (0,x,£). By choosing the test function to be Ao = fo(xixl + ^2^2)5 we have 

(Ao,EW £ )(0) = (/o,Ea £ + )(0). 
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Since (A ,EW e >(0) -> (A ,Wb)(0) = </o,Tr(n + W ))(0) and (f ,Ea%)(0) -> (/ ,Ea+)(0), 
we obtain Ea+(0) = Tr(IL).Wo(0)). The same discussion holds for Ea_ (0). 

In the first approximating martingale inequality, we choose s = 0, 

|(Ao,EWQ(t)-(Ao,EW e )(0)- f ds(($ + A)X ,EW £ )(s)\ < C Xo ,tV^- 

Jo 

By the choice of Ao, we further obtain 

<Ao,EW e )(t) - (A o ,EW £ )(0) - F ds{d t \ ,EW £ )(s) 

Jo 

=(fo,Ea%)(t) - (f ,Ea%)(0) - [ ds(d t f ,Ea%)(s) 

f-t 



^(/o,Ea+)(t)-(/ ,Ea + )(0)- / ds(d t f , Ea+)(s) 

Jo 

as e — > 0. 

For the term jj ds(>Uo, EW e )(s), we need to calculate x\A* {EW £ )x\ + X2A*(EW e )x2- 
After some lengthy algebra, we get 

xtA*(EW e )x! + x* 2 A*(EW £ )x 2 = - \ j- + +T{Ea £ + ,Ea e _) + (/), 

where (I) includes terms containing c?-,d?-. By the bound of cross modes, we check that 

Jq ds(fo, (J))(s) — )• as e — > 0. Thus, we have shown that (Ea+,Ea_) is the weak solution 
to 

d t Ea + H — — — = T(Ea + ,Ea_) 

with the initial condition given by Ea±(0) = Tr(n±Wo(0)). 

In the same way, if we choose Ao = foiuiUi + 2/2 2/2 ) 5 we can show (Ea + ,Ea_) also 
satisfies 

c£ • V x Ea_ 
<%Ea_ H — — — = T(Ea_,Ea + ). 

By the uniqueness of the solution to the above transport equations, the proof is complete. 

□ 



4.3 Convergence of the second moment 

We now prove the second approximating martingale inequality. Recalling the construction 
of jjb £ = fj,o + y/efJ,i >£ + e/^2,£, we need to prove that the correctors Ci, i = 1, . . . , 5 are small. 

For C4, we have the following lemma. 
Lemma 4.5. (C±,Ci){t) < C fl(h T£ d uniformly in t G [0,T] and AG V. 
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Proof. Recall that 

3 3 3 3 

L 2f i = ieE{- ^2 r°r fc ^ 1 i5b®Gi+F ®Gi r°T k fcsi-Yl f °f fe ^ 2 Fi®G +Fi(8)Go f °t k iC k e2 }. 

k=0 k=0 k=0 k=0 

Consider one term, for example 

JWx&fb ® Gi = £ (W I dP t% P) ^ P ' £ Fo{x,i- P -)\ ® G X . 
fc=o fc=o V ^ / 

Since 



Gi(i,s,^,0 



jf ^ L 7£)? e " rAl( '' P)ze (p)^o(C - f ) - G„(* + |)G t (p) 



e 



and Gt(p) = X]fc=o 7°7 fc ^ L fc(i>.P)> taking one term in the sum and ignoring constant, we 
have 



7 Rd (27r) d 2 y 

For any element in matrix E(J), we know that it is a linear combination of terms of the 
following form 

r°° r i(xi —x2)-p ~ p p 

(i)= dr dpe i R mn (r,p)f(x 1 ,£ 1 --)g(x 2 ,& + -)T 1 (r,&,p)T 2 (r,&,p), 
JO JR d z 1 

and /, g, T±,T 2 are from Fq, Go, e~ rAl ^ 2 ' p \ e~ r ^ 2 ^2.-p) respectively and we can assume they 
are all real. So 



r roc /*oo ^ ^ i( x l — X 2)*(P1 — V2) 

\\{i)\\r2(M2d\ = / / HR mn (s 1 ,p 1 )R mn (s 2 ,p 2 )e i dsids 2 dpidp 2 dxidx 2 d£id^ 
where 

# = n n ^2,6+^) n ^i.^**) n ^a^). 



By density argument, we can assume /(xi, £i)<7(x2, £2) — ^-1(^1 + ^2)^-2(^1 — ^2)^3(^1)^4(^2) 
for some good function hi, then we have 

h =\h l {x l +x 2 )\ 2 \h 2 { Xl -x 2 )\ 2 n ^(6-|) n Mk+f) n t ^ s i^2,pi) n ik^,*^. 

i=l,2 t=1.2 i=l,2 i=l,2 
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Change variables y\ = x\ + x 2 , y 2 = x\ — X2, and integrate in yi, we have 

[ h 3 (Ci-j) ] [ + ] [ ] [ ^( S2 ,6,p 2 ) 



L- 



c 



J 



t=l,2 



i=l,2 



i=l,2 



i=l,2 



-Rmn ( S 1 , Pi ) -R mn (s 2 ,p 2 )P( ^ ^ )dsid,S 2 dpi dp 2 d£ 1 d& 

for some constant G and z> is the Fourier transform of of 1 7^2 1 2 - Since Tj is bounded, 
integrating in £1 and £2 yields: 

f OO POO 



< 



J 



P2 — Pi 
I Rmn 

(si,p 1 )R mn (s 2 ,p 2 )i)( )\ds 1 ds 2 dp 1 dp 2 . 



Changing variable p 2 = p± + u, then integrating in pi,si,s 2 , since R mn £ S{M x R d ), we 



have 

Therefore we see that 

□ 



(«")!& 



|z>(— 



e . 



Since /ii = Fo <2> Gi + Fi <8) Go and /i2 = Fq ® G 2 + F 2 ® Go + /22, and we already have 
the control for Fi,Gi, F 2 ,G 2 , we only have to show 

(A2,e, M2,e) < (^Ho,T 

to conclude that (Gj, Gj) < G^^e for i = 1,2, 3, 5. The equation satisfied by /i2 is 

3 

' /i2 + imoc 2 (r + f°)/x 2 



fc=i 



SA2 - c 

3 

k=l 

L 2f i + «eE{^ r°r fc /C^F ® Gi - F Gi ^ r ^^ + ^ f °f fc /C^ 2 Fi ® G - Fi ® G ^ f°f A 



r°r fc p fel ^-^) + f°f fc p fc2 ^-^ 



im c 2 /i2(r + f ) 
3 



k=0 



k=0 



k=0 



k=0 



We just have to note that the matrices 
Ai :=!> 



k=l 
3 



r°r*jfc(« + + f°f*p fc * 2 (« + 1) 



im c 2 (r° + f ^ 



fc=i 



r°r fc p fcl ^-|) + f f fc p fe2 (ie-|) 



+ im c 2 (r° + f ^ 



are both of the form iQ for some real symmetric matrix Q, thus any element in the matrices 
e~ rAl , e~ rA2 is bounded. The rest of the proof is similar to the one for X 2>£ . 

With the second approximating martingale inequality, we can prove the following propo- 
sition. 
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Proposition 4.6. If(a^_,a £ _) converges weakly to (a+,a_) in C([0, T\; Hm ), then (a+,a_) 
is the unique weak solution to (2.9). 

Proof. First of all, we claim that a±(0,x,£) = Tr(II±Wb(0, x, £)). By the proof in Propo- 
sition 4.4, Eay-(0, sc,£) = Tr(IT-|-Wo(0, a;, £)), so we only have to show that a±(0) is deter- 
ministic, and this comes from the fact that 

(/o® 5o ,E{a + }®E{a + })(0) = (F ® G , W ® W )(0) = </ ® 5 o, ® a+})(0), 

since they are all the limit of (Fq ® Go, E{W £ ® W e }) (0) and we have chosen Fq = /of^ixf + 
X2X2), Gq = go(x±xl + X2X2). The same discussion holds for a_ if we choose Fq, Go with X{ 
replaced by yi. 

In the second approximating martingale inequality (3.6), let s = 0, 

\{F ® G ,E{W £ ® W e »(t) - (F G ,E{W £ ® W e })(0) 

- / cfe(F ® [(d t + A)G ] + [(d t + A)F ]®G ,E{W E ®W E })(s)\ < y/i. 
Jo 

and by the aforementioned Fq, Go, we have 

(F ® G ,E{W £ ® W £ }) = (/o ® gu,E{<4 ® <}) 

and 

<F ® [(ft + A)G ] + [(d t + A)F ] ® G ,E{W £ ® W £ }) 

= (/0 ® 5-0, E{Q+ ® (x*A*W e Xi + W e X 2 ) + (xlA*WeX! + 4i*f £ l2) ® «+}) 

+(/o ® d t g + d t fo ® goMa% ® a*.}) 
By the same discussion as in Proposition 4.4, let e — > 0, we have 
{fo®g ,E{a + ®a + })(t) - (/„ ® 50, E{o+ ® a+})(0) 
= / ^((ft + 4^f)/oC3 5o + /o® (ft + ^^f) 50 +,E{a + ®a + })( S ) 

+ / ds(f ® 5o,E{a + ® T(a+, a_) + T(a+, a_) £S> a+})(s). 

JO 

Note that we can define some operator such that a+ ® 7~(a+, a_) + 7~(a+, a_) ® a+ could 
be written as a functional of (a + ® a + ,a + ® a_,a_ ® a+). Therefore, we have derived 
an equation satisfied by (E{a + ® a + },E{a + ® q_},E{q_ ® a+}). By the result from 
Proposition 4.4, we check that (Ea + ® Ea + ,Fa + ® Ea_,Ea_ ® Ea + ) satisfies the same 
equation. 

By choosing other forms of Fo, Go, we can derive an equation system satisfied by (E{a + ® 
E{a + ®a„}, E{a_ ®a + }, E{a_ ®a_}). We check that the same system of equations is 
also satisfied by (Ea+®Ea_|_,Ea_|_®Eo:_,Ea_®Ea:_|_,Ea_(S>Ea!_), and since they share the 
same initial condition, the solution is unique, therefore we know (a+,a_) is deterministic 
and satisfies (2.19). The proof is complete. □ 
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4.4 Tightness 

In this section, we prove that (a £ + , a £ _) is tight in C([0, T]; Hm)- So together with Proposition 
4.4 and 4.6, we have finished the proof of the main theorem. 

Proposition 4.7. (a £ + ,a £ _) is tight in C([0, Tj; Hm)- 

Proof. By Proposition 2.5 we only have to show for any good function / that the process 
(/, afj_) G C([0,T];M) is tight. Take of, for example and note that 

(/,<) = (f(x 1 x* 1 + x 2 x* 2 ),W e ). 

Define Ao = f(x\x\ + X2X2) so we have (f,a £ + ) = (Ao, W E ). Recall that 

(A , W £ )(t) = f* ds((d t + A)Xq, W £ )(s) + Q]: £ AWe\{t) - (d,W e )(t) - f ds(C 2 + C 3 + C 4 , W e )(s) 

Jo Jo 
:= x £ (t) - y e (t). 

By our previous results, we have 

\y £ (t)\ = \(C 1 ,W £ )(t)+ f ds(C 2 + C 3 + C 4 ,W £ )( S )\<C ftT Vi 

Jo 

uniformly in t £ [0,T] and A £ V, so y £ (t) =^ in C([0,T];R). Then we know that if x £ {t) 
has a weakly convergent subsequence, so does {Xo,W e )(t), i.e., the relatively compactness 
of x £ (t) implies relatively compactness of (Ao, W £ )(t). Thus by the Prohorov theorem [8] we 
only have to prove that x £ (t) is tight in C([0,T];R). 

We apply the following sufficient conditions to show x £ (t) is tight. First, a Kolmogorov 
moment condition [8] in the form 

E{\x £ (t) - x e (ti)Hx e (ti) - x £ (sW} < C ftT \t - (4.2) 

with 0<s<ti<t<T and 7 > 0, /3 > 0. Second, we should have 

lim limsupP{ sup \x £ (t)\ > R} = 0. 

R^oo 0<t<T 

Since both {Xo,W e )(t) and y £ (t) are uniformly bounded, the second condition is auto- 
matically satisfied, so we only have to prove (4.2) for 

t 

l,er 





where ^' £ [lT £ ](t) is a martingale. We have 



x £ (t)= / d s ((d t + A)x ,w £ )( s ) + g L x ' o £ [w £ }(t), 



An 

E{\x £ (t)-x £ (s)\ 2 \F s }<2E{\ [ t dT((d t +A)Xo,W £ )(T)\ 2 \T s } + 2E{\g{ o £ [W £ ](t)-g{ o £ ^ 

J s 

■■= (I) + (II) 
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and (I) < C(t — s) 2 . To estimate (//), we only have to calculate the increasing process 
associated with the sub-martingale |^J, e [l^e](*)| 2 - We check 

\Q X x"[We){t)\ 2 = martingale part+- I ds (Q\(\ £ ,W £ }\ 2 - (X £ , W £ )(W £ , QX £ ) - (W £ ,X £ )(QX £ , 

£ Jo 

Moreover, the integrand Q\{X £ ,W £ )\ 2 - (X £ , W £ )(W £ , QX £ ) - (W £ , X £ ){QX £ , W £ ) is of order 
e by Lemma 4.2, 4.3 and the boundedness of Q. Therefore, we have (II) < C(t — s), which 
leads to 

E{\x £ (t) -x £ (s)\ 2 \F s } < C{t-s). 
In order to obtain (4.2), we note that 
E{\x £ (t) - x e (ti)H* B (ti) - x £ (s)P} = E{E{\x £ (t) - XeihWlFtJlxeih) - x E (s)P} 
<E{(E{\x £ (t) - x £ (h)\ Vt 1 }) i \x e (h) - x e (a)|T} < C(t - ii)iE{E{|:c £ (ii) - x £ (s)y\T s }} 
<C(t - £ 1 )iE{(E{|x e (<i) - x £ (s)\ 2 \T s }fy < C(t - - a)i < C(i - 

By choosing 7 € (1, 2), we finish the proof. □ 

5 Slow time fluctuations 

In this section, we discuss briefly the case when the temporal random fluctuations are slower 
than the spatial fluctuations of the potential. This leads to a kinetic regime with elastic 
collisions as the random fluctuations are now too slow to affect such a scattering; see [4] for 
a similar derivation in the setting of Schrodinger equations. 

More precisely, we consider the Dirac equation (2.9) with Ai{^, |) replaced by Ai(^, |) 
for some a £ (0, 1). The proof of convergence for a sufficiently large remains almost the 
same as before except that the infinitesimal generator of (W £ (t), A(t/e a )) becomes 

±E^ t {(X(V), W)}(t + h)\ h=0 = (±QX + d t X + A £ X, W). (5.1) 
So our test functions Ai,A2 in (3.10) and (3.12) become 

A i = " / e r£l - aQ e-^F ie - rA -drdp, 
J Rd (2ir) a J 

r iz-p /-co 

A 2 = " / / e " ~ aQ e-™{F 2 + £ £ X (27r) d 8(p))e- rA *drdp, 

where C £ Xo is defined by replacing R mn (r,q) in CXq by R mn (re l ~ a , q). Thus we see that 
Ai j£ is of order e a ~ l and A2, e is of order e 2a ~ 2 . To make those correctors C, small, we have 
to choose a G (f,l)- Note that for Ci, 62,64, a > ^ is enough, but for C3, it has to be 
greater than |. 

In the approximating martingale inequalities (3.4) and (3.6), we replace .AAo by the 
e— dependent operator | 2^fc=i( c 7°7 fe AfcAo + hDkXocy j k ) + C £ Xq. The rest of the proof is 



27 



the same except that when we pass to the limit, we need to calculate x\C*x\ + x\L* e X2 and 
this leads to some e— dependent scattering operator of the form 



r £ (a + ,a_)=-^ / (a+(q) - a+(fl) ^ u k (£, q )e a ~ 1 R kk ( cX+{q) 1 _ Q CA+(0 ,q - £)dq 
+ J f- E [ (a_(g) - a + (0) E g)^- 1 ^^ ^^^^ » g ~ Odq- 



Following the same type of proof in [4] and letting e — > 0, we arrive at the elastic scattering 
operator 

3 

T(a+) = e 2 / (a+(g)-a + (C))E^(^^( cA +(g)- cA +(O)^fc(0,g-e)^- (5-2) 

Therefore, in the slow temporal fluctuation case, when a G (f,l), we have the following 
transport equation system of the limit (a+,a_): 

<9*a± + C ^ = T(a±), (5.3) 

A ±(?J 

where T is the elastic scattering operator defined in (5.2). 

We see that the coupling between a+ and a_ appeared in (2.19) is inactive, and we 
expect (5.3) to hold in the limit of no time-dependent regularization, i.e., formally for 
a = 0. 

Remark 5.1. The condition that a > | could be relaxed somewhat. As a matter of fact, if 
we construct the test function as 

N 

En 
e 2 An )£ , 

n=l 

and follow the same procedure, we can show that A nj£ is of order e na ~ n , and l| + not — n > 
gives us a > |. For the corrector 

3 i 3 

o„,fc 



it is of the order e Na a - , and Na — > leads to 

1 1 

a> 2 + 27V 

Therefore, by expanding in higher order, we can relax the assumption to be a £ (^,1). 
Note that a = ^ corresponds to the regime considered in [10]. 

For slower time fluctuations of the media, i.e., when a < ^ or even a — > 0, other tech- 
niques than the Markovian regularization considered here presumably need to be developed, 
and the use of the diagrammatic techniques as in [11, 23] is currently unavoidable. 
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6 Conclusion and further discussions 



In this paper, we derived the kinetic limit of the Dirac Equation with time-dependent 
random electromagnetic field. We have shown that the cross modes cfj,dfj converged to 
zero weakly in space while the limiting propagating modes a + = an +022 and a_ = 611 + 622 
satisfied a transport equation system. In addition, the temporal regularization brings some 
new features (inelastic scattering) to the scattering structure, which disappear when the 
random fluctuations are slower in time. 

The method we use relies on the fact that the random field A(t,p) is Markovian in 
t. By constructing appropriate test functions, we prove some approximating martingale 
inequalities, and together with the tightness result, we pass to the limit. We should mention 
that our approach is restricted to the L 2 case in the sense that the initial data should be 
appropriately generated such that 11^(0) ||| 2 ( R d) is of order e d / 2 . In the setting of initial 
conditions ||^ ,£ (0)||^ 2 ( K d) of order O(l), we do not expect convergence of the energy density 
to a deterministic limit; see [6] in the setting of the scalar Schrodinger equation. 

Some generalizations could be obtained by the same approach. We could for instance 
show that the complete set of modes {af •, 6f •} converges weakly to the solution to a transport 
equation system of higher dimensions in (2.21). The case when R mn 7^ can also be handled 
by more strenuous computations. 

The same approach is expected to extend to linear hyperbolic systems with random 
coefficients such as those considered in [22]. This is currently under study. 
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